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FOREWORD 
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States  Air  Force  Contract  No.  AF  h9{638)-6k6,  Project 
No,  9751;  Task  No.  37510.  The  work  was  administered  by  the 
United  States  Air  Force  Office  of  Scientific  Research. 


SUMMARY 


The  stability  of  a  laminar  flame  subjected  to  a  small  perturba¬ 
tion  is  examined.  An  important  practical  case  is  considered  in  which 
the  amplitude  of  the  disturbance,  as  measured  by  the  flame  displace¬ 
ment,  is  small  in  comparison  with  the  wave  length  of  the  disturbance, 
but  not  necessarily  small  in  comparison  with  the  flame  thickness. 
Starting  from  the  basic  equations  governing  a  reactive  mixture  and 
making  use  of  an  expansion  procedure  in  the  parameter  X  ,  which  is 
defined  as  the  ratio  of  the  wave  length  to  the  flame  thickness,  the 
following  results  are  obtained: 

1,  For  £_  >>  1,  the  phenomenological  theory  of  Landau,  in  which 
the  flame  is  idealized  as  a  surface  of  discontinuity  propa¬ 
gating  with  a  constant  speed,  is  justified.  Accordingly, 
Landau's  conclusion  that  a  flame  is  unstable  to  small  dis¬ 
turbances  is  strictly  valid  for  disturbances  of  sufficiently 
long  wave  length, 

2,  The  second  term  in  the  expansion  in  l/oC  yields  a  correc¬ 
tion  to  Landau's  theory.  Both  the  coefficent  of  thermal 
conduction  and  viscosity  enter  into  the  correction.  The 
stabilizing  Influence  of  these  additional  terms  is  indicated 
the  details  are  reported  elsewhere. 


ON  THE  STABILITY  OF  LAMINAR  FLAME 


1.  Introduction  There  are  two  approaches  to  the  study  of  the  stability  of 
laminar  flame.  In  the  first  approach,  the  flame  is  treated  as  a  surface  of 
discontinuity  endowed  with  certain  properties  and  the  stability  of  flame  is  ex¬ 
amined  from  the  hydrodynamic  viewpoint  without  any  reference  to  the  detail 
chemical  reactions  involved.  This  is,  therefore,  a  phenomenological  approach. 

The  principle  contributions  of  this  school  are  represented  by  the  work  of 
Landau  (ref.  1)  and  Markstein  (ref.  2).  Landau  considers  the  flame  front  as  a 
surface  of  discontinuity  which  propagates  with  a  constant  speed,  small  in 
comparison  with  the  local  sound  speed,  and  across  which  there  is  a  discontinuous 
change  in  temperatm'e.  He  found  that  such  a  surface  is  unstable  to  distur¬ 
bances  of  all  wave  lengths.  Markstein  improved  Landau's  model  by  allowing  the 
flame  speed  to  depend  on  the  flame  curvature.  Markstein  found  that  the  flame 
is  stable  with  respect  to  disturbances  of  sufficiently  small  wave  length. 

The  second  approach  to  the  study  of  flame  stability  takes  into  account  the 
structiire  of  the  combustion  zone  and  is  certaiiiLy  the  more  correct  approach 
since  it  gives  a  detailed  analysis  of  the  changes  which  take  place  in  the  com¬ 
bustion  front.  The  principle  contributions  of  this  school  are  represented  by 
the  work  of  Rosen  (ref.  3)  Barenblatt  and  Zeldovich  (ref.  and  Menkes  (ref.  5)- 
However,  these  authors  have  only  examined  the  stability  of  flame  relative  to 
one  dimensional  disturbances. 

In  this  paper  we  shall  attempt  to  bridge  the  gap  between  the  two  approaches. 
If  the  stability  of  the  phenomenological  models  of  Landau  and  Markstein  are 
truly  relevant  to  the  study  of  flame  stability,  it  must  be  possible  to  justify 
these  models  from  the  first  principles.  In  other  words  we  must  base  our  analysis 
on  the  fluid  dynamic  equations  governing  a  reactive  mixture  and  consider  the 
problem  of  the  stability  of  flame  relative  to  three  dimensional  distui’bances . 
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It  is  found  that  Landau's  model  can  be  justified  provided  that  the  wave  length 
of  the  disturbance  is  long  enough,  in  compeirison  with  the  flame  thickness,  that 
the  changes  across  the  flame  are  Independent  of  the  detailed  profile  of  the 
flame.  Corrections  to  Landau's  model  for  disturbances  of  moderate  wave  length 
are  also  deduced.  It  is  shown  that  the  effect  of  heat  conduction  on  flame 
stability  can  indeed  be  described  by  a  model  similar  to  that  used  by  Marksteln. 
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Y  -  V.- 

~  N 


(2.6) 


(2.7) 


where  is  the  molecular  weight  of  the  ith  species,  since  the  molar  concen¬ 

tration  n^  is  related  to  the  mass  density  by 

(2.8) 

For  a  mixture  of  perfect  gas,  the  equation  of  state  is 

|o  =  vv-RT  (2.9) 


where  R  is  the  universal  gas  constant  and  n  is  the  total  numbers  of  moles 

M 

of  all  the  species  in  the  mixture,  i.e.,  n  =  ^  . 

'■  = » 

The  motion  and  changes  of  state  of  a  chemically  active  mixture  are  gov¬ 
erned  by  the  following  equations : 

i)  Continuity  Equation 

ii)  Mass -Balance  Equation  for  Eac.i  Species 

iii)  Equation  of  Motion 

iv)  Energy  Equation 

v)  Equations  of  State 

The  derivation  of  these  equations  can  be  found  elsewhere  (e.g.,  ref.  6  or  7). 
The  resvilts  are  summarized  below 
Continuity  Equation: 


=  0 


(2.10) 


where  u  is  the  mass -average  velocity  of  the  mixture,  t  denotes  time  and 
V  is  the  del  operator. 
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Mass  Balance  Equation  for  Each  Species 


= 


(2.11) 


where  ^  denotes  the  total  derivative,  W;_  denotes  the  total  rate  of  production 
of  the  ith  species  per  unit  volume  and  is  the  diffusion  velocity  of  the 
1th  species.  The  diffusion  velocities  are  related  to  the  binary  diffusion 
coefficients  Dlj  by  the  relation: 


(2.12) 


i  i  \ ,  i  , .  M 


Here  we  have  ignored  the  pressure  and  thermal  diffusion.  Only  N-1  of  the 
above  equations  are  independent,  since  the  mass -average  velocity  u  of  the 
mixture  is  defined  such  that 


2-  X  'i  =  o 


(2.13) 


The  conservation  of  mass  demands  that 


^  W-  =  O 


Equation  of  Motion 


9  =  '7-  b 

^  ut.  r 

where  the  stress  tensor  p  has  the  typical  component 


(2.14) 


(2.15) 


f,.,  --  -  (.Sj,  *  ^  .  I  5,^1^.]  (2.16) 

In  (2.l6),  ^  is  the  coefficient  of  viscosity  of  the  mlxtiure,  Is  the 

unit  tensor  and  summation  convention  has  been  employed 


-4- 


Energy  Equation 


(2.17) 

where  H  is  the  enthalpy  of  the  mixture  given  by  (2.5),  ^  is  the  heat  flux 

vector  and  ^  is  the  viscous  dissipation: 

Ip  =  (hi*'  ■ 

(2.18) 

The  heat  flux  vector  Q  consists  of  two  parts:  that  due  to  conduction  and 
that  due  to  diffusion, 

a,  .  -\'7T  +  4  (£.19) 

where  \  is  the  coefficient  of  thermal  conduction  and  is  the  enthalpy  of  the 
ith  species. 


Equation  of  State 

In  terms  of  the  mass-fraction, 

N 


the  equation  of  state  is 


(2.20) 


3.  Chemical  Reactions  Involving  Two  Species  Obviously  the  above  system  of 
equations  is  very  complicated.  To  simplify  somewhat  the  analysis,  we  shall  limit 
our  discxission  to  the  special  case  of  a  single  reaction  involving  two  species: 

M,  and  M^.  We  assume  that  they  enter  a  reaction  of  the  form 


(3.1) 


where  k  ,  k  denote  respectively  the  forward  and  backward  specific  reaction 
f  b 

rate  constants.  Clearly,  -  50|ni,  =  ^^.l)  includes  two  types  of 

reactions  of  practical  importance: 

i)  Unimolecular  reversible  reaction  in  which  case 
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11)  A  one-stop  decomposition  reaction  In  which  case 


k  -  O 

b 


Y-  - 


In  our  case,  N  =  2  and  It  Is  convenient  to  Introduce  a  new  variable  defined  by 

\  =  Xc  +  "'4.  ^ 

where  Is  the  Initial  (unburned)  mass-fraction  of  the  species  It  follows 

from  (2.1)  and  the  relation  - that 

V  V/  ,1-  ll.xb) 

where  Yo  Is  the  Initial  mass  fraction  of  the  species  M,.  Is,  therefore. 


a  measiare  of  the  composition  of  the  mixture.  We  shall  denote  the  rate  of  reaction 
(3.1)  by  fiTj  which  Is  directly  proportional  to  the  rate  of  Increase  of  the 
mole  number  per  unit  mass  per  unit  time.  The  precise  definition  of  F  Is  embodied 
In  the  defining  formula 


Clearly,  by  (2.14),  we  have  also  the  similar  formula 


(3.3a) 


(3.3b) 


The  Law  of  Mass  Action  states  that  If  (3.1)  represents  an  elementary  reaction 


f  =  (3.4a) 

where  and  are  known  functions  of  temperature.  If  terms  of  5 


and  T  we  have 


kjd)  r  9CY.o^>n,yjX)]'^'_  t  j] 

s  J  ^  ^  J  (3.1 


In  a  unlmolecular  reversible  reaction  we  have; 


f  =  k^tT)  [■?!?.( ^  i<(T))\ 

where  KCx)  =  equilibrium  constant. 


(3.4c) 
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When  there  are  only  two  species,  (2.12)  and  (2.12)  can  be  combined  to  give 


the  simple  relation 

-J) 


(3.5) 


where  5)  *  3^ij.  =  Di,,  ■  Introducing  (3.5)>  (3.3)>  (3.2)  into  the  Mass -Balance 
Equation  for  the  species  (i.e.,  ei^.  2.11),  we  obtain  the  simple  equation 


S  M  -  =  fF 


(3.6) 


We  next  consider  the  energy  equation  for  the  case  Let  us  substitute 

(2.5)  and  (2.19)  into  (2.17).  Introducing  the  specific  heat  at  constant  pressure 
Cpj  of  the  species  and  making  use  of  (2.11)  we  obtain  an  alternate  form  of 
The  Energy  Equation 

Z?v,Ceal?-=;vTj- If  (3.T) 

i  -  \ 

Let  us  now  introduce  the  simplifying  assumption  that  Cp;  can  be  approximated 
by  a  mean  Cp  which  is  the  same  for  both  species  over  the  range  of  temperature 
involved.  This  Implies  that 


H,  =  Cp  Ct-To)  +  Hco 

E<1.  (3.7)  can  then  be  reduced  to  the  simple  form; 

fCp5I  .  |>  =  7.(A.vT)+$  4  fAF 


(3.8) 


(3.9) 


where 


(3.10) 


Finally,  by  substituting  (3.2)  in  (2.10),  the  equation  of  state  assumes 
the  simple  form 


(3.11) 
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where 


=  ■R  £ 

•  mi 

is  the  gas  constant  of  the  unburned  mixture. 


(3.12) 


Summarizing,  in  the  case  of  a  single  reaction  involving  two  species,  there 
are  five  variables  which  describe  the  state  of  the  fluid:  Ty  w.  a»ie.l  ^ 

These  variables  are  governed  by  the  five  equations:  (2.10),  (3.6),  (2.15),  (3.9) 
and  (3.11).  With  a  set  of  properly  imposed  boundary  conditions,  the  five 
equations  determine  completely  the  five  unknowns. 


4.  Laminar  Flame  We  consider  a  low  speed  steady  one-dimensional  flow  with  a 
rapid  change  of  temperature  occuring  in  a  narrow  zone  (width  L)  in  which  most 
of  the  chemical  reactions  occur.  By  low  speed  it  is  meant  that  the  Mach  number 
is  small  in  comparison  to  unity.  Such  a  flow  field  describes  the  change  of  state 
in  a  plane  lamineir  flame.  Once  we  have  chosen  the  solution  for  such  a  flow,  we 
will  be  interested  in  its  stability.  We  will  assume  a  small  perturbation  of 
the  one- dimensional  flow  and  see  if  there  is  amplification  or  decay.  The 
problem  will  be  supposed  two-dimensional,  but  the  perturbation  being  slight, 
the  governing  equations  may  be  linearized. 


5.  One-Dimensional  Steady  Flow  The  equations  governing  such  a  flow  are: 


r  u  ^  -I  K 


^0*^0  “  Coni 


.■fani" 


(5.1) 


(5.2) 


(5.i^) 
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In  order  to  proceed  further  we  must  specify  the  chemical  reaction  taking  place 
in  the  flow.  For  our  purpose,  it  is  as  general  and  more  convenient  to  specify 
an  adequate  temperature  distribution  and  then  deduce  the  production  term 
from  (5.3).  Several  such  temperature  distributions  might  be  convenient.  Thus, 
we  may  have 


"  Tooc 


o  +  e 


(5.5) 


This  is  shown  in  Fig.  1.  Or,  we  might  also  take  a  linear  distribution  of 
temperature  Inside  the  flame,  although  such  a  distribution  is  certainly  un¬ 
realistic  at  the  limits  of  the  flame.  (See  Fig.  2)  From  experimental  data,  it 
can  be  concluded  that  the  Reynolds  nimiber  of  the  flow  based  on  the  "flame  width" 
is  of  the  order  of  unity.  This  can  also  be  seen  from  (5.3)  if  one  assumes  that 

oCoUn  ^nd  are  of  the  same  order,  which  is  certainly  the  general 

JO  p  o  '■dx^ 

case.  If  L  is  the  thickness  of  the  flame,  then  I  ~  and 


\ 

U 


That  is,  l^eTr  ~l  where 


H  ooo  Op 
^000 


(5.6a) 

(5.6b) 


Now  for  an  ordinary  gas,  Tr~l  so  that 


IPe  ~  I  (5.6) 

Eq.  (2.1)  then  show  that  ~ 

1=  =  t=.-  OCM-jJ  (5.6) 

Since  the  Mach  number  of  a  flame  is  small  in  comparison  with  unity,  it  follows 
from  (5.6)  that  the  pressure  is  constant  throughout  the  flow.  Consequently, 
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the  pressure  and  viscous  terms  in  (2.3)  can  be  neglected  in  comparison  with  the 
conduction  term.  We  have  them 


(5.7) 


which  is  of  the  same  form  as  (2.i4-)  with  T  replaced  by  g  If  the  Lewis  nvmber, 


Le 


^  P 
A  ' 


(5.8) 


is  chosen  to  be  equal  to  one^ it  follows  from  (2.4)  and  (5.7)  that  §  =  T  + 

a  constant,  vhlch  can  be  determined  by  the  condition  at  infinity.  Thus,  as  ^  ^ 
Tj, Tot-5  and  |  O  .  Hence 


In  particular^  if  the  distribution  (5.5)  is  valid 


(5.9) 


5  ■ 


From  the  equation  state  (assuming  a  uni-molecular  reaction), 


(5.10) 


-t- 

«0*O 

u 


(5.11) 


so  that 


i^o  -  3'.«L'+  4  (5.11) 

for  the  distribution  (5.5).  The  velocity  distribution  is  given  by  (5.2); 

l^c  =  To  (5.12) 

*0« 

Finally  is  given  by  (5.4)  or  by  (5.7) 


6.  Some  General  Remarks  on  The  Perturbed  Flow  The  perturbed  flow  is  ob¬ 
viously  governed  by  linear  equations.  Since  the  unperturbed  field  is  one-dimensional. 
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the  coefficients  in  these  equations  will  depend  only  on  x  and  not  on  ^  ^  and 

t.  The  stability  of  flame  can  either  be  treated  as  an  initial-value  problem 
or,  more  simply,  by  a  formal  analysis  in  which  one  harmonic  component  is  con¬ 
sidered.  Without  any  loss  of  generality  we  can  consider  the  solution  as  in¬ 
dependent  of  ^  .  Furthermore,  since  the  coefficients  of  the  equations  are  In¬ 
dependent  of  y  and  t.  we  may  assume,  in  a  formal  aneilysis,  a  dependence  in  y 
and  t,  of  the  form  where  W  is  the  wave  number.  In  other  words,  all 

the  perturbed  quantities  will  be  of  the  form  -(ca)  ^ 


Let  us  consider  the  points  of  the  flame  at  a  certain  temperature.  It  can 
be  easily  seen  that  they  will  be  distributed  on  a  line 

Censtfvot-  ^  (6.1) 

where  i)  and  the  constant  are  functions  of  temperature.  If  ^  denote  the  actual 
temperature;  T.  the  undisturbed  temperature  and  T  the  perturbed  temperature, 

•5^  *  To  +  T  (6.2) 

Let  be  a  point  on  the  constant-temperature  surface  T  .  At  a  later  in¬ 

stant,  it  will  be  located  at  the  point  X.  where 

C(.y)  =  To  IXo)  (6.3) 


For  small  disturbances  5:~Vo  and  we  have 

'Qcxi  ^  ^CXo)+  U-^o-)pciCo)  (6.^) 

Solving  for  'X  and  making  use  of  (6.2)  and  (6.3),  'we  have 

X  =  (6.5) 
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_  0  vl<W-t-LO“t 

But  I  CXo’)  -  jCA-o)  e  and  from  (6.3)  we  can  express  as  a  function 

of  ^  .  Hence,  (6.5)  can  be  written  as 

X-  +  <^(^)  (6.6) 

which  proves  the  statement  (6.1) 

The  at  a  certain  reference  temperatvire  represents  the  amplitude  of  the 
disturbance.  For  linearization  to  be  possible,  we  must  have 

1>  <<  k~‘  (6.7) 

From  the  kinetic  theory,  it  is  known  that 

-V,  (6.8) 

where  X  is  the  mean  free  path  and  a,  the  sound  speed.  Now  I  ,  therefore 

L  -  4  (6-9) 

But  our  equations  are  valid  only  if  we  deal  with  a  continuous  flow;  that  is,  all 
distances  entering  the  problem  must  be  much  greater  than  i  ,  Eq.  (7»9)  shows 
that  L »  ,  if  M  <«•  I  ,  which  is  our  case.  But  it  also  shows  that  no  distances 

can  be  allowed  to  be  of  an  order  smaller  than  L  .  in  particular,  D  is  of  order 
L  or  greater,  i.e.,  the  flame  displacement  is  at  least  of  the  same  order  as 
the  thickness  of  flame.  Consequently,  it  is  not  possible  to  linearize  the 
equations  inside  the  flame  using  X  as  a  variable.  O'.e  way  to  overcome  this 
difficulty  is  to  use  the  temperature  as  the  Independent  variable  instead  of  X.  , 
This  is  particular  suitable  since  the  undistrubed  field  is  also  given  in  terms 
of  the  temperature. 
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7.  Transformation  of  Variables 


We  now  introduce 


The  limits  between  the  zones  are  not  well  defined,  but  from  (5.5) 


cIt 


—  _ 


eo  L 


t-T-TooJ  L'^^-Tcoc  C‘'*‘  ^00^3 

so  that  we  cein  also  say: 
upstream: 


flame: 


V  downstream: 


4^0 

<1T 


tiT 


T -  Te  JO  Tec 

L 

Tsoo 


»-  - _  »  Lx 

T-T.^C>  +  'C«ec'i  T^eo 


(7.6) 


(7.7) 


Although  this  definition  is  equivalent  to  the  previous  one  if  the  distribution 
(5.5)  is  valid,  it  is  neveT'if.gUss  preferable,  since  physically  the  temper¬ 
ature  distribution  must  be  continuous,  but  the  condition  of  continuity  is  less 

cjy 

stringent  upon  the  derivative  ^  which  can  vary  fast  from  a  near  infinite 
value  outside  the  flame  to  a  finite  value  inside. 


8.  Continuity  Equation  The  continuity  equation. 


4-  (^*1) 

it  dx  ^^4 

expressed  in  terms  of  the  new  variables  s 


it  or  at  "aT 


dT 


(8.2) 


and  using  the  linearization  process: 


clx„  _  'djc 

at;  3t  di  it 


+  +  4 


O  8X0 

J”  dr 


0 


(8.3) 


9.  Momentum  Equation:  The  two  momentum  equations  are 


fill 


&x 


+  ^r'^]  -  - 


+i 

rt  3  ix*-  3 


0*  i 


1 


(9.1) 
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^x 


L  3 


j. 

3.  'ixb'i 


J)X''  J 


In  terms  of  the  new  variables,  we  have 

- = -If; 

s  I  f?  I?  -  i?  H  *  4^  *  j  =  - If  T 


(9.2) 


where  "V.T."  stands  for  the  vlscoixs  term.  These  terms  are  calculated  below; 

,"^1^  ^  =  'is  -  lif- 
^.av  ox*-  ■»!  St'-  ST  ST*- 
i 


Jt  lx  Villi'  _  lx  ^  ^  "Sx 

'aT'  “  '-»T^  aTa>i  ot'-Vt  a>j  at  sT&Tbi^  Sr'-'Sj 

_  "ix  iii— (Is  "ii  H 'Is ‘L£  lii  bSt 
LSt/  ^3*'’  *  a^(  iTai|^'-»T*  ar'-iT  VTSi^ 

-V  ..  .  1 


(9.3) 


etc. 


yty.'-a  /'Si'x  \  1  pxN^a'X  1  _  "ix  liX/^ 

St  'aTx  ■'■  S'^ViT/ata^  '■H 


^x) 


Applying  the  linearization  procedure,  we  obtain  finally  the  perturbed  momentum 
equations : 


Ijff  [f “-k  *  #i  '  h  li  '  “.1^  '  ^'k)  *  HI 

*  ‘If  art' 

■  ^4•Ax,3■^^<  _lt!<ix,  du,  ^  du,  Mx.\  4.  1  /Ay.S  Sv  1 

’’’  p  dr  aT‘~  dT  3t*-~  3t  a'V'’ ^  a>j>^  ~  d"^*  I"! '■cTf/  IvaiVara'^J 

,  is^f  f  a  I  M  er  \  4  It  dxo  -  d|-e  is  1 

I  'fj/  ''  at-  ht  ^  ^4T  j  0|  r  jr  J 

=  U  r  .  'll.  ■Ls  -  +  -  I^K.A^  4  _l  Jt,..  3y; 

r  L  3  ''dT/  ax  ax^  cV'f’  ai  3>^A^J  «toii  j  cjt*-  33  ir>? 


»y^ 


-  1  aV  T 

<'t  <At  J 


(9.4) 


where  use  has  been  made  of  the  relation 


U.  -  ^ 


OOo  ^“ 

'o 

\«»» 


10.  Diffusion  Equation  The  diffusion  equation  for  a  two  component  mixture 
with  Lewis  number  equal  to  unity  is 


^  [U  x,,3  ^X'-M  1  -  v-^7r)  =  i>T 

J  L  at  ape  ati  J  (  '*' 

a  .15. 


(10.1) 


In  terms  of  the  new  variables  we  have 


(10.2) 


»  1  ^  fv ^  /  ^- 1  +  A.  r  ^  1  ll  \  !lS  _  ^  '  111 

CpaTL^iT/  irj  ^  Cf  L  at  atj^  ?>tJ  L  arJ 

T  Cf,T,  g)  is  a  known  function  of  j’ ,  T  Q.%d  ^ 

When  the  flow  field  is  perturbed  slightly 

F  -  Pc  i  f; 

F.  -  f(|fj/  Ksil 

where  =  Y  §„■)  'the  production  x-ate  at  T  in  the  undisturbed  flow. 
Remembering  that  *aTo  linearized  diffusion  equation  assumes  the  form: 

f42?ol'  r  4B--  )  -  p  -ir  P  F  "ifo  ?)  P  F 

VdY^  'dT  dr  ‘<5?  5r if  ar/  -FocTt  n  ~  ^ 


(10.3) 


dT 

X  2.^  dx,  p  u  cl^. 
^  ^  d]  ^  dr 


(10.4) 


_  "a'*?  4^',  tllr  dxo  5I  1  51  -  '^^“1 

Cf  -  a-T’dT  "  aT*<Fr  “  <57^ 5T  '3?;  b>i>  '■Yi/  ■»■«»  JtoJ 


(11.1) 


(11.2) 


T  ^6 

11.  Energy  Equation:  The  Energy  Equation: 

^  -t  ^  d  pAf 

expressed  in  terms  of  the  new  variable  is 

tl-xf  1  J4  - 

-\  r  .  ax  y\  ax  ,■^p^x 
~  L  oT*'  OT  ara-^  a>^ 

niff  Li  *  pflFj 

Now,  the  dissipation  function  ^  is  given  by 

S  /  it  r  (^f  +  ^  -»•  f 

2/r-  =  5L^axJ  +  U^;J  35^^ 

Thus  it  contains  terms  of  the  following  types: 

W  ,  ,  5u  ^  ■  L  ^  ■air  .  ( 'bir- 

ax  dy  I  ax/  t  a^/  -*  v  ax-'  ax  j  '  ^ 

Obviously  in  the  linearisation  process  only  the  first  two  terms  will  give  any 


(11.3) 


contribution. 

Thus  we  may  write 


where  is  the  steady  dissipation  function,  and  ‘t*,  the  first  order  dissipation 
function: 

(11.5) 

^  3  3-?aTI1-  "  3  s-pVa-T-J  J 

v. 


(ii.J^) 


A  ij  f-i  lr+1  di(.4?«'3u  +  iax  1 

<13,  ■  t  a  ar  3>7  3  3t  <*T  3  67  Var  <  J 


Applying  the  linearization  procedure  to  (11.2),  we  obtain  the  energy  equation: 


(11.6) 


12.  Independence  Principle  and  Justification  of  Landau's  Model  Since 
the  order  of  magnitude  of  each  physical  quantity  can  be  easily  estimated,  it  is 
possible  to  rewrite  our  equations  such  that  each  term  in  its  nondimens ional 
form  will  be  of  the  order  of  unity.  It  is  convenient  to  introduce  the  following 
parameter: 

£  is  a  measure  of  the  intensity  of  the  distrubance  since  it  is  the  ratio  of 
flame  displacement  to  the  wave  length  of  the  distrubance.  In  oirr  problem,  we 
have  three  length  scales:  k"'  ,  D  and  L.  In  the  -direction  we  must  ob¬ 
viously  use  W  *  as  the  nondimens ionalizing  factor;  in  the  x-direction  we 
should  make  a  distinction  between  the  perturbed  and  unperturbed  quantities: 
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we  must  use  L  to  non-dlmenslonalize  x  and  D  to  non-dimens ionalize  x.  In 

O 

order  to  facilitate  the  identification  of  the  terms  we  will  keep  the  same 
notation  whether  the  quantity  has  a  dimension  or  not.  Table  I  will  be  use¬ 
ful  in  establishing  the  relation  between  the  dimensional  and  non-dimensional 
form  of  the  same  variable.  When  the  governing  equations  are  written  in  the 
non-dimensional  form,  the  following  dimensionless  numbers  appear  in  the 


Variable  Quantities 

U. 

■U. ,  IT  ,  CO 

1 

X 

So 

s 

% 

h 

fo 

$ 

~T 


Table  I 

To  Be  Non-dimens ionallzed  By 

£.^00. 

L 

k-' 

CpTop.  A'' 
eCpTopo  A"' 

Mobo/ 

g  Uc..  / 


equations:  Mach  number,  Prandtl  number,  Reynolds  number,  (Lewis  number,  which 
has  been  taken  equal  to  unity).  In  addition,  we  have  the  dimensionless  number 
X=tkL)-'  (12.2.) 
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I 


which  is  the  ratio  of  flame  thickness  to  the  wave  length  of  the  disturbance,  and 
the  dimensionless  nxomber  (/V  =  which  represents  the  ratio  of  heat 

of  reaction  to  the  heat  capacity  of  the  mixture.  In  the  following,  we  shall 
neglect  terms  of  the  order  in  comparison  with  unity.  In  terms  of  the  non¬ 
dimens  ionalized  variables,  simimarized  in  Table  I,  the  governing  equations 
assme  the  following  form: 


JT '■  8't>'  ^  4^  T-x  zyy  ^ 


aV_  .  j./42fofiV  _  I  ^  1 

^ 


(12.3) 


(12.4) 


'•af/  L  7t  ^  STt  j:  +  3^  J 

_J—  f  it  X  ^  4  X  iHa.  .  J-  ^  dtf] 

i?€  XHdfj  n’aj  -^sTJiidT  5tj  (12.5) 

[5  -  f-s-  M M ‘  ‘  fei?- 

■  ■  ■flr'il'  ■  fer?  (12.6) 


(12.7) 


J  = 


(12.8) 


Eq.  (12.3)  is  the  continuity  equation;  (12.4)  and  (12.5)  are  the  momentum 


equations;  (12.6)  and  (12.7)  respectively  are  the  diffusion  equation  and  the 
energy  equation  ;  (12.8)  is  the  equation  of  state,  all  written  in  dimension¬ 
less  form.  Note  that  (12.8)  shows  that  the  compressibility  effect  is  negligible. 
The  non-dimensional  variables  f’.  ,  5*0  ,  w-o  ,  f ,  ^  describing  the  undis¬ 

turbed  flame  are  related  to  one  another  by 


dx 


=  =  si  (•-!.)  = 

Jx  dxvfo^ 


_  i£  J-  duo  dx, 

m'ldxJ  -  <ri-  ^dxj 


(12.9) 


(12.10) 


Substracting  (12.7)  from  (12.6)  and  making  xise  of  (12.9),  we  obtain: 


(^o)‘r  4X0 +11  1  5.  -1-  r  (12.11) 

Lj;.  2't  dT  ax  -I  fftPr  •- iT»-JT  dx'ax  dl  -1 


Now  the  solution  of  this  equation  having  an  exponential  dependence  on  ')^  and 
'C  which  satisfies  the  homogeneous  boundary  conditions  C  *  o  at  the  upstream 
and  downstream  infinity  is  clearly 

?  =  O  (12.12) 

This  equation  states  the  familiar  fact  that,  when  Lewis  number  is  equal  to  one, 
x.he  composition  of  the  mixtiure  is  directly  related  to  the  temperature  of  the 
mixture. 


It  is  evident  that  a  flame  can  be  approximated  as  a  surfa.ce  of  discontinuity 
only  when  the  wave  length  of  the  disturbance  is  much  greater  than  the  thickness 
of  the  flame.  Hence,  Landau's  model  is  strictly  applicable  only  when  X . 

In  Landau's  model,  the  flame  is  considered  as  a  surface  of  discontinuity  across 
which  there  is  a  discontinuous  temperature  change  and  which  propagates  with  a 
constant  speed  relative  to  the  fluid  immediately  ahead  of  the  surface.  If 
Landau's  phenomenological  model  gives  a  correct  description  of  a  flame  when 
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X  -?  cc  j  it  must  be  possible  to  justify  it  by  using  the  fluid  dyneimic 
equations  governing  a  reactive  mixture  as  the  basis  of  analysis.  Furthermore, 
if  Landau's  model  can  be  justified,  it  paves  the  way  for  Improving  the  model 
so  that  the  effect  of  heat  conduction  and  viscosity  may  be  taken  into  account. 

To  justify  Landau's  model,  let  us  expand  the  various  dependent  variables 
in  an  Inverse  power  series  of  the  parameter  X  • 


(12.13) 


etc.  Substitute  these  expressions  into  (12.3),  (12.4),  (12.5),  (12.12),  (12.7), 
(12.8)  and  equate  terms  of  the  same  degree  in  X  . 

First,  let  us  consider  terms  of  order  X'.  Eqs.  (12.4)  and  (12,7)  yield 
the  following  relations: 


X  ■= 

X.  +  ^  ^  1 

1=  = 

F.  i  K 

U  X 

Cl  0  4  U, 

-  <i>!o  -t  =  O 

ai  jiT*-  dT«-  op  ^ 

'  4  O  '^0  I  1  -  —  P  o  W  1  s  Q 


(12.14) 


These  equations  can  be  considered  as  a  pair  of  simvLLtaneous  homogeneous  algebraic 
equations  governing  the  two  variables:  and  .  Since  the  deter- 

mlnent  of  the  system  is  equEil  to  -  (gy^^  F,  fo''^  1®  different  from  zero 

inside  the  combustion  zone,  it  follows  that 

Tt  =  °  (12.15) 

Consequently,  when  the  wave  length  of  the  disturbance  is  large  in  comparison 
with  the  thickness  of  the  flame  so  that  all  terras  of  order  and  higher  can  be 
neglected,  the  thickness  of  the  flame  remains  constant  when  it  is  slightly  dis- 
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turbed.  In  other  words,  the  effect  of  the  disturbance  is  to  carry  the  flame 
along  without  changing  its  structure. 


Let  us  next  collect  terms  of  order  X* 


From  (12.4),  we  find 


df.  ft;  1.^)  ^  =  O 


(12,16) 


Using  (12.15),  we  may  rewrite  (12.16)  as 


-  c, 


(12.17) 


where  C,  is  a  constant  which  may,  hop^ever,  be  a  function  of  and  T  .  Eq. 
(12.17)  expresses  the  fact  that  the  mass  flux  through  the  flame  is  constant. 
From  (12.6),  we  find: 

(fcfe  [  V.  *  •  4  [-  ih  hi -*  "’Vp  - 


Now,  from  (12.9),  I  •  This  equation  can  thus  be  rewritten  as 


which  may  be  integrated  at  once  to  give 


(12.19) 


[7.-*  1^-J  >  C. 


(12.20) 


The  equation  has  a  si’^file  meaning.  First  of  all,  we  note  that  is  the 

slope  of  the  constant  temperature  surface  so  that  represents 

the  tangential  velocity  of  the  flow  relative  to  an  observer  fixed  on  the  con- 
stant-temperatxrre  surface.  The  term  1  thus 

represents  the  velocity  gradient  at  the  constant -temperature  surface  and 
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eq.  (12.20)  states  that  tangential  momentum  transport  minus  the  shearing  stress 
is  continuous  across  the  flame.  If  the  viscous  effect  at  two  sides  of  the 
combustion  zone  is  negligible,  we  have  simply  the  condition  of  continuity  of 
tangential  velocity: 


V  Wo 


(12.21) 


We  next  turn  to  eq.  (12.7).  Collecting  terms  of  order  ^  'w®  have 


This  equation  can  be  interpreted  as  £in  equation  governing  the  change  in  the 
thickness  of  flame.  For,  if  eq.  (12. I7)  is  substituted  into  (12.22),  we  obtain 
a  second  order  differential  equation  for  X,  : 


_L_ 

5T>- 


3  ^  ^^<0  ] 

<»T  J 


(12.23) 


Obviously  the  change  in  the  flame  thickness  caused  by  distxjrbances  in  the  flow 
must  depend  on  the  reaction  mechanism  and  the  structure  of  the  fame,  as  is 
indeed  embodied  in  the  terms  and  Xo  in  (12.23).  This  effect  will  be  im¬ 
portant  only  when  JL  becomes  less  than  unity. 

Eq.  (12.22)  can  also  be  regarded  as  an  equation  governing  the  change  of 
(u.,  -  ^  )  through  the  flame.  For  I  7  expect  that  the  variation 

of  all  quantities  of  the  order  of  X  or  even  of  the  order  X.  should  be  in¬ 
dependent  of  the  structure  XqCt)  of  the  flame.  Assuming  this  "Independence 
Principle",  we  conclude  immediately  from  (12.23)  that 


dT 


•=  o 


and  hence 


^0 

5't 


(12.24) 


(12.25) 


-23- 


Equation  (12.25)  states  that  the  constant -temperatxire  surface  moves  at  a 
constant  velocity,  since  the  perturbed  velocity  relative  to  an  observer  fixed 
on  the  constant-temperature  surface  is  u^-  which,  according  to  (12.25),  is 
zero.  This  substantiates  Landau's  model. 

Eq.  (12.24)  states  that,  even  to  the  order  of  ,  the  flame  thickness 

is  independent  of  the  flow  disturbances. 

If  we  had  applied  the  Independence  Principle  to  eq.  (12.20)  or  (12.19), 
we  would  have  obtained  directly  (12.21). 

Finally,  let  us  turn  to  eq.  (12.4).  Collecting  terms  of  the  order  of 
and  making  use  of  (12.15),  (12,24),  (12.25),  we  obtain; 

^  (12.26) 

That  is,  the  perturbed  pressure  is  continuous  across  the  combustion  zone. 

Obviously  eqs.  (12.25),  (12.17),  (12.21),  (12.26)  are  identical  to  the 
four  boundai’y  conditions  which  can  be  deduced  for  a  surface  of  discontinuity 
of  Landau's  type.  This  justifies  Landau's  model  for  a  flame  when  the 
Independence  Principle  is  applicable.  Of  course,  this  principle  itself  is  a 
postulate,  the  correctness  of  which  must  be  verified  by  the  internal  con¬ 
sistency  and  by  deductions.  However,  one  thing  is  clear:  when  this  principle 
is  not  applicable,  Landau's  model  is  inadequate  to  describe  a  flame  front. 
Moreover,  it  will  be  seen  that  Important  deductions  can  be  made  from  this 
hypothesis  which  defines  mathematically  meaningful  problem  and  provides  physically 
meaningful  results. 

13.  Effect  of  Conduction  on  Flame  Stability  To  deduce  a  correction  to  Landau's 
model  for  not-so-large  X.  1  ^6  shall  compute  the  terms  of  order  X  in  the 
expansion  (12.13).  We  shall  first  consider  the  effect  of  heat  conduction  on  the 
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stabilityj  and  we  shall  neglect  the  effect  of  viscosity  altogether.  Collecting 
terms  of  the  order  of  '/[£,  and  applying  the  Independence  Principle,  the  con¬ 
tinuity,  momentim  and  energy  equations  yield  the  following  system  of  equations: 


(13.1) 

(13.2) 

(13.3) 
(13. J+) 


To  be  exact,  eq.  (13.1)  is  derived  from  (12.3);  (13.2)  from  (12.4);  (13.3)  from 
(12.5)  and  (13.4)  from  (12.7).  In  addition,  we  have  (12.24): 

^  Q  (13.5) 

Eq,  (13.2)  and  (I3.I)  may  be  combined  to  give 

-f  =  O  (13.6) 

aT  0-r 

where  use  has  been  made  of  (13,5)  and  (12.9).  This  equation  may  be  used  in¬ 
stead  of  (13,2).  Now  (see  eq.  12. 15),  it  follows  that  eq.  (13.4)  im¬ 

plies  (13.1).  Consequently,  eq.  (I3.I)  is  not  an  independent  relation  and 
need  not  be  considered  further.  Finally  eq.  (13.3)  can  evidently  be  written  as 


(13.7) 


Eqs.  (13.4),  (13.5),  (13.6),  (13.7)  govern  the  changes  of  u, ,  V,  ,  ,  x, 

in  the  flame.  Now  ^,,^0  K  '  ^0  satisfy  eqs.  (12,25),  (12.21),  (12.26) 
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and  (12.15).  Thus,  they  also  satisfy  the  following  equations: 


H  ° 

is  4  ^  0 

ar 


If  ve  multiply  (13.4),  (13.5),  (13.6),  (13.7)  by 
equations  of  (13.8),  we  obtain: 


'b'x 


It  4  =  O 

*T  ST 


^  J  ^  O 


ST 


(13.8) 


and  add  to  the  corresponding 


(13.9) 

(13.10) 

(13.11) 

(13.12) 


which  reduce  to  the  boundary  conditions  associated  with  Landau's  model  if  we 
let  X  “>  00  .  Eq.  (13.10)  shows  that  X.  is  independent  of  T  so  that  the 
flame  thickness  is  not  affecteJby  the  disturbances  in  the  flow.  Eq.  (13.9)  shows 
that  the  flame  speed  is  affected  by  the  curvature  of  the  flame.  The  negative 
sign  in  (13.9)  corresponds  to  the  fact  that  the  flame  speed  increases  when  the 
flame  front  is  concave  toward  the  unburned  gas .  Finally,  (I3.II)  and  (13.12) 
can  be  integrated  to  give  the  boundary  conditions: 


[3'  4  au  = 


au." 


V  -  V 

''t  -  Vt 


(13.13) 

(13.14) 


where  the  prime  and  double  prime  signify  the  state  of  flow  at  two  sides  of  the 
flame  and  Is  the  tangential  velocity 


v;  -  v  4- 

T 


(13.15) 
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The  boundeiry  conditions  (13.9)  -  (13.12)  are  identical  to  those  obtained  by  a 
different  method  in  an  earlier  analysis  by  the  authors,  (ref.  8)  There  it  is 
shovn  that  these  boundary  conditions  are  applicable  under  a  less  stringent 
condition.  It  is  further  shown  that  the  conduction  effect  stabilizes  the  flame 
if 

I  _L 

URei;  ^  (13.16) 


14.  Effect  of  Viscosity  on  Flame  Stability  We  next  examine  the  effect  of 
viscosity  by  the  same  method.  Here  we  shall  ignore  the  influence  of  heat  con¬ 
duction.  Collecting  terms  of  the  order  of  Vji  and  applying  the  Independence 
Principle,  the  continuity,  momentum  and  energy  equations  yield  the  following 
system  of  equations. 


d. 

dT 


Lfo(«.  - 


O 


IF, 

bT 


L?.  N  =  r  J.  3^^ 

J  1?e  L 


lx.  Juo  _ 
■ST  a-x 


^*1 


o 


1  -1- 

R’e  bTO)j 


1 


In  addition,  we  have 


IF.  = 


O 


(li^.l) 

(14.2) 

(li^.S) 

(14.4) 


(1‘^.5) 


Since  (l4.4)  and  (l4.5)  imply  the  relation  (14.1),  the  latter  need  not  be  con¬ 
sidered  fvurther.  Substituting  (l4.4)  eind  (12.21)  into  (l4.2),  we  obtain 


_  4  j_  ^ 

bT  "  5  -We.  2^*- 


(14.6) 
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Finally,  substituting  (12.25)  and  (12.15)  into  (1^.3)  we  have 


(14.7) 


■where  use  has  been  made  of  (l4.5).  Eqs.  (14-.4)  -  (14.7)  govern  the  changes  of 

combustion  zone.  No'w  i^o  /  t’o  /  <-  satisfy 

eqs.  (12.25),  (12.21),  (12.26)  and  (12.15)  so  that 


tt 


a-t 


(14.8) 


-2l  r 

4 


Vo+  u 


IS, 
o  a 


r![“J  =  c) 


Let  us  multiply  (l4.42)  to  (l4.7)  by  X  and  add  them  to  the  corresponding 
equations  of  (l4.8) .  We  obtain: 
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These  boundeiry  conditions  reduce  to  that  associated  with  Landau's  model  if  we 
let  X  60  .  They  have  been  applied  in  ref.  8  to  study  the  effect  of  vis¬ 
cosity  on  the  stability  of  flame.  It  is  shown  that  viscous  effect  cannot  statil; 
the  fleime  if  fias  a  vodwa  greater  ■than  one  .  Thus,  conduction  has  a 

greater  influence  than  the  viscous  effect  in  stabilizing  the  flame. 
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